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Abstract. We have made a new calculation of the cosmic 
ray electron spectrum using an anomalous diffusion model 
to describe the propagation of electrons in the Galaxy. The 
parameters defining the anomalous diffusion have been re- 
cently determined from the study of nuclei propagation. The 
predicted electron spectrum is in a good agreement with the 
measurements. The source spectral index, found from exper- 
imental data, in this approach turns out to be equal to 2.95. 



1 Introduction 

Observations of non-thermal radiation of the Galaxy stimu- 
lated investigations of propagation of cosmic ray electrons 
through the interstellar medium. Since basic paper of Sy- 
rovatskii (1959), the problem of calculation of electron spec- 
trum was considered in series of papers (see, for example, 
Ginzburg and Syrovatskii (1964); Berkey and Shen (1969); 
Shen (1970); Bulanov et al. (1972); Cowsik and Lee (1979); 
Nishimura et al. (1979); Berezinski et al. (1990); Atoyan et 
al. (1995); Moskalenko and Strong (1998)). The normal dif- 
fusion equation for concentration of the electrons with en- 
ergy E, 

N(r, t, E), generated by sources distribution with density 
function S(r,t,E), 

BN 8 

— = DAN(r, t,E) + — (b(E)N(r, t, E)) 



S(r,t,E), (1) 



has been used to study the electron energy spectrum modifi- 
cations in the interstellar medium (ISM). In the equation (1) 
D is the diffusivity, b(E) — the energy-loss rate of electrons. 
Recently, in the papers (Lagutin et al. (2001a,b); Lagutin 
and Uchaikin (2001)), new view of the cosmic ray propa- 
gation problem was presented. It has been shown that the 
"knee" in the primary cosmic ray spectrum is due to large 



free paths ("Levy flights ") of cosmic rays particles between 
magnetic domains — traps of the returned type. As the "Levy 
flights " distributed accordingly to inverse power law oc Ar~ 3 ~ a , 
r — > oo, a < 2, is an intrinsic property of fractal structures, 
in the fractal-like medium the normal diffusion equation (1) 
certainly does not hold. 

Based on these arguments in (Lagutin et al., 2001c) an 
anomalous diffusion (superdiffusion) model for describing 
of electrons transport in the fractal-like ISM was proposed. 
This superdiffusion equation for concentration of the elec- 
trons without convection has been presented in the form 



ON 
~dt 



-D{E,a){-A) a/2 N(r,t,E) 



dE 



(b(E)N(r,t,E))+S(r,t,E), (2) 



Correspondence to: Strelnikov (sdv@theory.dcn-asu.ru) 



where D(E, a) is the anomalous diffusivity and (—A)"' 2 is 
the fractional Laplacian (called "Riss" operator, Samko et al. 
(1987)). 

The solution of superdiffusion equation (2) in the case of 
point impulse source with inverse power spectrum and the 
behaviour of energy spectrum of electrons in high energy re- 
gion were found. 

The main goal of this paper is to calculate the spectrum 
of electrons from sub-Ge V to TeV energies in the framework 
of anomalous diffusion model . We don't use the assumption 
that the mean time of particle staying in a trap is finite. In this 
paper, similarly to (Lagutin and Uchaikin, 2001), we suppose 
that a particle can spend a long time in a trap, that is q(t) oc 
Bt~P~ l , t — > oo, (3 < 1 ("Levy trapping time "). The 
spectrum of electrons both in ISM and in the solar system is 
found too. 



2 Flux of electrons from point source 

The flux of electrons, J(r, t, E), is related to the source 
S(r , t , E ) by the propagator G(r, t, E; r , to): 



oo t 

J(r,t,E) = ^-J droj dE f 



dtn 



x G(r, t, E; r , t )S(r o ,t , E )S(t - 1 - r). (3) 



Here 



E„ 



dE 1 



(4) 



5(t — to — t) reflects the law of energy conservation in the 
continuous losses approach. 

The propagator in the anomalous diffusion model under 
consideration has the form (see, Lagutin and Uchaikin (2001)) 



G(r,t,E;r 0i to) = (D(E,a,0)tP) 3/a 

x¥t P) (\r\(D(E,a,^y 1/a 



(5) 



where 



tf < a ^> (r) = / q { 3 a) (rrl 3 ) q[ P - l) (r) r 3 ^ a dr. (6) 



Here q^ (r) is the density of three-dimensional spherically- 
symmetrical stable distribution with characteristic exponent 
a < 2 (Zolotarev et al. (1999)) and qf'^ (t) is one-sided 
stable distribution with characteristic exponent (3 (Uchaikin 
and Zolotarev , 1999). The parameters a, (3 are determined 
by the fractional structure of ISM and by trapping mecha- 
nism correspondingly, the anomalous diffusivity D(E, a, (3) 
— by the constants A and B in the asymptotic behaviour 
for "Levy flights" (A) and "Levy waiting time" (B) distribu- 
tions: 

D(E, a, j3) ex A(E, a)/B(E, (3). 

The energy-loss rate of relativistic electrons is described 
by the equation (see Atoyan et al. (1995)) 



dE 
~dt 



b(E) = bo + b 1 E + b 2 E 2 



b 2 {E + E 1 ){E + E 2 ), (7) 



where 6 = 3.06-10 _16 -n (GeV-s^ 1 ) is for ionization losses 
of electrons in ISM with number density n (cm~ 3 ), b^E with 
bi = 10~ 15 • n (s _1 ) corresponds to the bremsstrahlung en- 
ergy losses, and b 2 E 2 with b 2 = 1.38 • 10~ 16 {GeV ■ s)- 1 — 
synchrotron and inverse Compton losses, for B w 5/iG and 
w w l(eV/cm 3 ), E x w b /h, E 2 w b Y jb 2 . Using (7), the 



solution of the equaion (4) relative to E can be presented in 
the form (Atoyan et al., 1995) 



Eo(r) = 



E + E x 



l-{l-e- b ^){E + E 2 )/{E 2 -E l ) 



-El 



Taking into account that b\T < 3.15 -10 8 • n ■ 10 5 = 3.15 ■ 
n ■ 10 -3 Cl,we derive from the later equation 



Eo(t) 



E + Ei 



l-b 1 T(E + E 2 )/(E 2 -E 1 ) 



-El (8) 



With help of the equations (3, 5, 8) it's easy to calculate the 
flux of electrons for the sources interesting for astrophysics. 
For example, for point impulse source 

S(r, t, E) = S E- p 6(r)Q(T - t)@(t), 



Q(x) = 



1, x > 0, 
0, x<0, 



we have 



J T (r,t,E) = — So 

Ait 



min(t,l /b 2 (E+E 2 )) 

J dTE (r)- p 

max[0, t-T] 

x (X(E, t^)t^ 1 )- 3 /«(1 - b 2 r(E + E 2 y 2 

x^ { t (i) ^{KE^y- 1 )- 1 '*) , (9) 



where 



E (r) 



X(E,t)= J 



D{E\aJ) 



It should be noted that in the case (3—1, the equation 
(9) comes to the solution, obtained earlier by (Lagutin et al., 
2001c). If a = 2, (3 — 1, we have the standard solution. 

3 Energy spectrum of electrons 

The flux J of electrons due to all sources of Galaxy can be 
separated into two components: 

lkpc R 



lkpc 

J L (r < lkpc) + J G (r > lkpc). (10) 

The similar separation is frequently used in the studies of 
cosmic rays (see, for example, Atoyan et al. (1995) and ref- 
erences therein). 

The first component (L) in (10) describes the contribution 
of the nearby sources (at distance r < lkpc) to observed 
flux J. The second component (G) is the contribution of the 
distant sources (r > lkpc) to J. 
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Fig. 1. Energy spectrum of electrons near the solar system. 



The nearby sources used in our calculations is presented 
in Table 1. Based on result (9) we suppose 



Jl= J2 Mri,U,E), 



(11) 



(r<lkpc) 



where injection time T w 10 4 -r- 10 5 y. 

The second component (G) is evaluated under assumption 
that the distant sources (r > lkpc) are distributed uniformly 
both in space and time in the Galaxy. 

The parameters defining the anomalous diffusivity and used 
in our calculations have been recently derived from the study 
of nuclei propagation (Lagutin et al., 2001b): a — 1.7; j3 — 
0.8; D(E,a,/3) = D (E/lGeV) s with D w (1 H- 4) • 
10~ 3 pc 1 y~°' s and 6 = 0.27. Only one parameter p defin- 
ing injection spectrum of electrons in the sources is found by 
fit. Extensive calculations show that the best fit of experi- 
mental data may be get at p « 2.95. 

The spectra of L- and G- components and the total spec- 
trum in ISM are demonstrated in Fig. 1 . 

To describe the influence of the solar modulation , the 
force model of Axford and Gleeson (1968) is used: 



Jmod(E) 



E 2 - (m e c 2 ) 2 



[E + $(t)] 2 - {m e c 2 } 



;J[E + *(t)], 



where <!>(£) = 600MeV. The modulated spectra are shown 
in Fig. 1 too. 



4 Conclusion 

We have made a new calculation of the cosmic ray electron 
spectrum using an anomalous diffusion model to describe the 
propagation of electrons in the Galaxy. The parameters defin- 
ing the anomalous diffusion have been recently determined 
from the study of nuclei propagation. The predicted electron 
spectrum is in a good agreement with the measurements. The 
source spectral index, found from experimental data, in this 
approach turns out to be equal to 2.95. The proximity of this 
exponent to one obtained earlier (Lagutin et al., 2001b) for 
nuclei components (p « 2.9) can indicate the same mecha- 
nism of particle acceleration. 

We have shown that the sources of high-energy electrons 
(E > 100 GeV), observed in the solar system are relatively 
young local sources (r < 200 pc, t ~ 10 5 y), injecting 
particles during the time T ~ 10 4 -j- 10 5 y. The behaviour 
of spectrum in the low-energy region is defined by distant 
(r > lkpc) sources. 



Table 1. List of the nearest SNR 



Name 



r, pc 



t, 10 J y 



Source 



Lopus Loop 
Monoceros 
Vela 

Cyg. Loop 
CTB 13 
S 149 
STB 72 
CTB 1 
HB21 
HB9 
Monogem 
Geminga 
Loop I 
Loop II 
Loop III 
Loop IV 
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